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We analyze the ground state of the two-component gas of trapped ultracold fermionic atoms. 
We neglect the forces between atoms in the same hyperfine state (the same component). For the 
case when the forces between distinguishable atoms (i.e., atoms in different hyperfine states) are 
repulsive (positive mutual scattering length), we find the existence of critical interaction strength 
above which one atomic fraction expels the other from the center of the trap. When atoms from 
different components attract each other (negative mutual scattering length) the ground state of the 
system dramatically changes its structure for strong enough attraction - the Cooper pairs built of 
atoms in different hyperfine states appear. 

PACS Number(s):05.30.Fk, 03.75.Ss 
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I. INTRODUCTION 

After the first experimental achievement of quantum 
degeneracy in potassium ^"K gas ^ followed by the cool- 
ing of other elements (^Li) below the Fermi temperature 
U , the interest in ultracold atomic Fermi systems has in- 
creased significantly. Since at very low temperatures the 
s-wave scattering in a single-component spin-polarized 
Fermi gas is excluded by the statistics, the experimen- 
tal realization of quantum degeneracy requires actually 
trapping of two kinds of atoms. In fact, in Ref. the 
simultaneous trapping of different spin states of ^"^K is 
reported as a way which enables cooling the fermions 
by rethermalizing collisions undergoing between atoms 
in different spin states. Another idea is to employ the 
sympathetic cooling technique in a mixture of fermionic 
and bosonic isotopes like ^Li and ''Li or even in two- 
species Fermi-Bose mixtures like ^Li - ^'^Na Q or '*°K - 
**^Rb il]. 

Several aspects of the physics of degenerate Fermi gases 
have been investigated in the past few years. In Ref. 
the properties of a harmonically trapped spin-polarized 
(i.e., non-interacting) Fermi gas is studied by using the 
Thomas-Fermi approximation. This work was extended 
in several ways. One possibility is to consider the dipole- 
dipole forces as in Ref. j^, where the ground state of 
fermionic dipoles in the normal phase is analyzed and 
mechanical instabilities for large enough number of par- 
ticles or the dipole moment are discovered. Another way 
is going beyond the semiclassical approximation by in- 
cluding the discrete nature of the trap levels as well as 
by introducing the second component and the interaction 
between atoms from different components 

Not only static but also dynamic properties of cold 
Fermi gases have been investigated. For example, the 
collective excitations in degenerate fermionic gas are dis- 
cussed within the hydrodynamic approximation in Ref. 
and based on the sum rules in Ref. 9J. There has 
been developed an idea of generating solitons and vor- 
tices in a one-component Fermi gas in a normal phase 
by using the phase imprinting technique The pres- 
ence of vorticity in a trapped gas rotating at low angular 



velocity is analyzed in Ref. [lTI |. 

Much emphasis has been also put on the achievement 
of superfluid phase transition. There was a proposal 
to reach the Cooper-paired state at temperatures com- 
parable to the Fermi energy via a Feshbach resonance 
[l2| for a short-range interaction between distinguish- 
able fermions. It is plausible that the fermionic atoms 
confined in an optical lattice could also undergo the 
phase transition at high temperatures (one tenth of the 
free-space Fermi energy) independently of the scattering 
length |l3j |. The possibility of p-wave pairing through 
the Feshbach resonance for a short-range interaction 
in a single-component Fermi gas was discussed in Q 
whereas the achievement of the superfluid transition in 
dipolar gases (again in a one-component case) in Ref. 



Structure of binary Bose-Einstein condensates have 
been studied extensively over the last years Illlia. In 
Ref. |0| possible classes of solutions within Thomas- 
Fermi approximation are given and verified by numeri- 
cal integration of the coupled Gross-Pitaevskii equations. 
For the clarity only the case when all atoms repel each 
other is considered. It was found that for some range 
of parameters (scattering lengths) the ground state of bi- 
nary system is a symmetry-breaking solution. Moreover, 
it turns necessary to go beyond the Thomas-Fermi ap- 
proximation to obtain asymmetric ground state because 
the contribution of the kinetic energy (neglected within 
Thomas-Fermi approximation) is substantial in this case. 

The similar approach but for a two-component Fermi 
gas is reported in Ref. where the role of s- and 
p-wave interactions is investigated within the Thomas- 
Fermi approximation. The existence of the critical parti- 
cle number was shown for pure repulsive s-wave interac- 
tion above which the spatial separation of both compo- 
nents appears. However, after the separation is reached, 
the Thomas-Fermi approach does not yield the unique 
answer. Further assumptions are necessary. For example, 
consideration of configurations with minimal interface as 
energetically favorable leads to asymmetric ground state 
of two-component Fermi gas. For attractive s-wave in- 
teraction Ref. predicts the collapse of the system for 
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large enough number of atoms similarly to the case of 
Bose-Einstein condensate. 

In this paper we analyze the ground state of two- 
component trapped ultracold Fermi gas and demonstrate 
the existence of several characteristic coupling strengths. 
To this end, we use the mean-field Hartree-Fock method 
and also BCS related approach when necessary (in the 
case when different kinds of atoms attract each other). 
So, opposite to Ref. [TH, we go beyond the Thomas- 
Fermi approximation. On the other hand, we are re- 
stricted by numerics to rather small sample of atoms. 
Moreover, in this paper we present results of only one- 
dimensional calculations. 

The paper is organized as follows. In Sec. ^ based 
on the Lagrange formalism, we derive the equations 
which govern the two-component gas of trapped ultra- 
cold fermionic atoms. Sec. IIIII defines the mixed two- 
particle density matrix and calculates densities of each 
component. In Sees. II VI and IVl we discus the properties 
of the ground state of the system within the Hartree- 
Fock and BCS approaches. Finally, we conclude in Sec. 
EH 



II. EQUATIONS DESCRIBING 
TWO-COMPONENT FERMI GAS AT ZERO 
TEMPERATURE 



The time-dependent many-body Schrodinger equation 
can be derived as the Euler-Lagrange equation from the 
following Lagrangian density 



2 dt 2 dt 2m 



N 



considered as a function of the many-body wave func- 
tion ^'(ri, Tat) and its spatial and time derivatives. 
At low temperatures both components can be treated as 
non-interacting Fermi gases since the s-wave scattering 
is absent for spin-polarized fermions. The only interac- 
tion left is the repulsion or attraction between atoms of 
different spins. At zero temperature the wave function of 
two-component {N -f TV atoms) Fermi system is assumed 
to be the product of two Slater determinants 



^'(Xi, ...,XAr;yi, ...,yAr) 



(2) 



Inserting the wave function |5J into the Lagrangian ob- 
tained based on the Lagrangian density and integrat- 
ing each term over the appropriate {N — 1) coordinates 
for each component one obtains the effective Lagrangian 
density 
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which depends on one-particle orbitals of both types and 
their spatial and time derivatives. Deriving Eq. 
we have assumed the contact interaction between distin- 
guishable fermions with g being the coupling constant. 
Like in experiment T] both atomic fractions move in 
harmonic potentials and the ratio of frequencies is equal 
to oj^^-'/w^^^ = -^/O/T. Time-independent version of the 
Euler-Lagrange equations originating from the effective 
Lagrangian density J^J is given by 



fc2 N 

N 

El (1)|2 (2) (2) (2) 

\Vj \ - 4 'Pi 

t^l,2,...,N (4) 
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with ef''^ and e|^^ being the single-particle Hartree-Fock 
energies. The total energy of the system is expressed in 
the form 

(1) , J2) 



<vI/|H|*>=^ef)+^£, 
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'Pi (XAf) 



f N i^N) 



/|v.r^(x)f|^f(y) 



dx dy . 



(5) 



We solve numerically the set of Eqs. ^ to find the 
ground state of two-component degenerate fermionic 
gases. 
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III. ONE- AND TWO-PARTICLE DENSITY 
MATRICES 

Mixed two-particle density function is defined as 
pW(^)(xi,yi;x{,y{) = 

< ^Wt(x,)v3(2)t(y^)v}(i)(x;)V3(2)(y;) > (6) 
or equivalently in the position representation 
p«(^)(xi,yi;xi',y0 = 
y"*(xi,X2, ...,XAr;yi,y2, ...,yAr) 

x**(x{,x2,...,xjv;y{,y2, -^Yn) 
X(ix2...dxjv (iy2-..dyjv (7) 

One-particle densities are derived from mixed two- 
particle density functions 1617(1 by contracting over x or 
y variables 

p(^)(xi;xO=<^«nxi)^(^nxO> 

pi'^(yi;yi')=<V'''^^(yi)^^'nyi)> (8) 

or equivalently 

Pi (xi;xi) ^ ^(xi,yi;xi,yi) dyi 

P?\yi;yi) = / 4'^^'^(xi,yi;xi,y0 dxi (9) 

In the Hartree-Fock state given by |5J) two-particle 
density matrix, since single-particle orbitals are orthog- 
onal, reduces to 

^E^r^*(xi)^r*(yi)^«(x0^r(yO (lO) 

and is always a product of one-particle density matrices 
(two-particle correlation function vanishes) 

(1)/' '^ (1)*/ ^ (1)/ '^ 

Pi (xi;xi) ^j;^2^'Pl i^im (xi) 

i 

p'iHyi-M)-^f:^?'*{yi)^r\yi) (n) 

The density of each component equals the diagonal part 
of corresponding one-particle density matrix and for the 
Hartree-Fock state (|2J is given by 

1 ^ 

/^''^(-i) = ^Ei^r'(xi)p 

i 

'^^'^(yi) = ^Ei^f'(yi)i' (12) 

i 

In Section we show that for strong enough attraction 
between atoms better description of the ground state is 
BCS like one with significant correlations developed. 



IV. POSITIVE SCATTERING LENGTH 

In this Section we consider the ground state of one- 
dimensional two-component Fermi system with repulsive 
forces between atoms of different components. We solve 
numerically the set of Eqs. Q (in its one-dimensional 
version) by using iterative procedure. At zero step we 
pick up the Hartree-Fock orbitals, for example as the 
eigenstates of the pure harmonic potentials for each com- 
ponent - these functions are called zero-step orbitals. 
Then we calculate potentials which come from the pres- 
ence of other component and build Hamiltonian matri- 
ces for both parts of the system. Matrix elements are 
calculated in sets of basis functions which include corre- 
sponding harmonic oscillator eigenstates. The number of 
basis functions is typically much bigger than the number 
of particles A''. At each step both matrices are diago- 
nalized and new single-particle wave functions are ob- 
tained as well as single-particle energies. This procedure 
is continued until all quantities remain unchanged with 
prescribed accuracy. 
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FIG. 1: Energy of 3 4- 3 atoms system as a function of pos- 
itive interaction strength obtained witliin the Hartree-Fock 
approach. For weak enough repulsion (g < 4) all initial con- 
ditions lead to the same density profile which is the density 
of the ground state. For strong repulsion the density pro- 
file of the ground state changes qualitatively its shape (see 
Fig. |5J. Various sets of zero-step orbitals give different local 
minima of energy. Insets show the zero-step densities used 
in the numerics: (a) gives the results marked by solid line, 
(b) corresponds to dotted-dashed line, and (c) leads to both 
dotted and dashed lines. See text for further description of 
the procedure. 

We summarize our results in Fig. ^ There exists the 
critical value of the interaction strength (5 w 4 for 3-1-3 
atoms system) above which the density profile of the sys- 
tem changes qualitatively its shape (see Fig. [SJl. It is 
profitable from energetic point of view that the compo- 
nent confined in a lower frequency trap is repelled from 
the center of the trap. Although the potential energy is 
increased in such a way, the total energy, however, is de- 
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FIG. 2: Density profiles of tiie ground state of two-component 
Fermi gas for various positive interaction strengths. Left and 
right columns show the case of 3 + 3 and 4 + 4 atoms systems, 
respectively. For small coupling constant {g < 4) both frac- 
tions almost fully overlap each other whereas beyond the short 
intermediate interval (around g = 4) fractions separate. Full 
curves correspond to the first (higher frequency trap) com- 
ponent. The broken curves indicate the second component 
confined in a lower frequency trap. 



creased because the interaction (repulsion) is weakened. 
The solid line in Fig. ^ shows the energy of the ground 
state of two-component system. Its density is presented 
in Fig. [3 (left column). The fact that the density is 
not symmetric with respect to the center of the trap (the 
lowest frame) is explained by the presence of odd number 
of atoms in each fraction. For example, for 4 -f 4 atoms 
system the density of the ground state is symmetric (see 
the right column of frames in Fig. O . The interplay be- 
tween the potential and the interaction energies of the 
gases is even more subtle than just described. It turns 
out that the symmetry of the ground state strongly de- 
pends on the geometry of the traps as shown in Fig. |3 
for a particular value of the coupling strength. The verti- 
cal dotted line indicates the ratio of the trap frequencies 
used throughout this paper. Depending on this ratio the 
ground state of the system switches between the asym- 
metric density profile (boxes in Fig. |3J) like in the middle 
frame in Fig. 2]and symmetric one (circles in Fig. O) 
as in the upper frame in Fig. 01 The vertical dashed 
line marks the geometry (w^^-'/^^^''— 0-945) in which the 
ground state of the system becomes degenerate. We have 



checked that increasing the repulsion strength between 
the components the ratio w'-^-'/'^''^'' that favors the degen- 
eracy of the ground state is getting lower. Eventually, for 
g > AA the ground state of 3 -I- 3 atoms system becomes 
asymmetric under the considered trapping potentials. 

For large enough g we find several local energy minima 
depending on the chosen set of zero-step orbitals. The 
lowest minimum is accessible provided one uses for the 
second component harmonic oscillator eigenstates which 
are first cut at the center and then pulled out as a zero- 
step orbitals (inset (a) in Fig. If plain harmonic os- 
cillator eigenstates are taken as the initial wave functions 
(inset (c) in Fig. ^ one follows the dashed line and even- 
tually ends in a domain structure state (the lowest frame 
in Fig. 21). The dotted line matches the dashed one in 
Fig. Hand was obtained by solving the Euler-Lagrange 
equations corresponding to the effective Lagrangian den- 
sity (121 under the assumption of adiabatic approxima- 
tion. Finally, the dotted-dashed line in Fig. ^represents 
the local minimum with the densities of both fractions be- 
ing shifted with respect to each other (for large enough g 
- the middle frame in Fig. 0J) . This minimum is obtained 
when at zero step of iteration one takes shifted left (for 
one component) and right (for another component) har- 
monic oscillator eigenstates (inset (b) in Fig. ^ . Nothing 
changes qualitatively when one goes to larger systems as 
shown in Fig. [S] which is the case of 50 -I- 50 atoms sys- 
tem. The dense structure visible in the lowest frame of 
Fig. [Slindicates that for much larger number of atoms we 
should expect a smooth overlapping distribution of both 
components in the corresponding case. 
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FIG. 3: The energy of an asymmetric density profile states 
(boxes) and symmetric ones (circles) as a function of the ratio 
cj'-^'/a;'-^' for 3-1-3 atoms system. The vertical dotted line 
shows the geometry used in this paper. 



V. NEGATIVE SCATTERING LENGTH 

In this Section we consider the ground state of two- 
component system with attractive forces between atoms 
of different components. First of all, we have checked 



5 




FIG. 4: Density profiles for tlie coupling constant g = 10. 
The upper frame corresponds to the ground state whereas 
the middle and lower ones to the excited states which show 
qualitatively different patterns; separation (the middle one) 
and the domain structure (the lowest one). 



that in this case the Hartree-Fock energies and densi- 
ties of the ground state do not depend on the choice of 
zero-step orbitals. The resulting ground state" energy 
is plotted as a dashed line in Fig. and the densities 
of both fractions (for g = ^10) are shown in the inset 
(also by dashed line). It is clear that for strong enough 
attraction between atoms of different fractions a new den- 
sity pattern with visible fragmentation is developed. Fur- 
ther analysis reveals the appearance of energy gap in the 
single-particle Hartree-Fock spectrum (see Fig. [7|l . The 
energy gap separates off the populated Hartree-Fock lev- 
els. If one of the atoms occupying the highest level is pro- 
moted to the state just above the energy gap the total 
energy of the system is increased and can be calculated 
based on formula (Q. It turns out that the surplus in 
energy is of the order of the energy gap. Existence of 
such an energy gap in a single-particle spectrum could 
be a signature of appearance of pairing phenomenon in a 
system. This possibility will be verified in the following 
subsection. At the same time no such gap is developed 
when distinguishable atoms repel each other (Fig. |S1 the 
upper frame). Note, however, that an energy gap is also 
developed (Fig. |S1 the lower frame) in the case of repul- 



FIG. 5: The sequence of frames similar to that in Fig. |l]but 
for the system of 50 -I- 50 atoms and g = 37. Solid curves 
show the densities of the first component (loaded into higher 
frequency trap) whereas the dashed ones correspond to the 
second component (loaded into lower frequency trap). The 
lowest frame is, in fact, the case of very dense domain struc- 
ture. 



sion between different atoms in antiferromagnetic phase 
present in excited states ( as in Fig. 0] the lowest frame) . 



A. BCS formulation 

To check the idea whether the energy gap developed 
in a single-particle energy spectrum is connected with 
the appearance of pairing phenomenon in a system, we 
have reformulated our problem. Within Hartree-Fock 
approach each particle interacts with all the others in the 
second component. Now, we allow for formation of pairs 
including atoms of different spins. More precisely, only 
configurations in which different atoms populate single- 
particle states of the same quantum number n are per- 
mitted. Hence, the ground state wave function is written 
in the following form [T^ 

\'^BCs>=Y[{un+Vnatbt)\0>, (13) 

n 
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INTERACTION STRENGTH (osc. units) 

FIG. 6: Ground state energy of 3 + 3 atoms system as a func- 
tion of interaction strength for both the Hartree-Fock and 
BCS approaches. For strong enough attraction the ground 
state of the system shows the appearance of Cooper pairs. 
The inset shows the density (in oscillatory units) calculated 
via the Hartree-Fock (dashed line) and BCS (solid line) meth- 
ods. Note that the densities for both fractions look the same. 



where ui 



1 for each n (this assures the normal- 



ization condition) and operators ci^,b^ create fermions, 
respectively, of first (a+) and second (5+) component. Of 
course, is the probability that the pair state ipiPipi? 
is occupied. 

The Hamiltonian of the system includes kinetic and 
potential energies of both parts as well as the interaction 
between distinguishable atoms and is given by 



71, k 



H = 

+ Y.<m\El^l + v}Z\l>b+h 

m,/ 

+ ^ < nm\Vint \ kl > 



(14) 



n,m,k,l 



The mean value of the Hamiltonian in the ground state 
given by ((T^ is calculated as 

EbCS - E(^" ' + ^" ^)^" + E ^nn^nnvl 



n, m 
n ^ m 



(15) 
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< nn\Vint\ mm > 
— < nm\Vint \ nm > 



(16) 



The first term in (|15|l is the sum of kinetic and poten- 
tial energies of all atoms, the second and third (the first 
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FIG. 7: Single-particle Hartree-Fock energies for various 
strength of attraction between atoms of different components. 
For strong enough attraction the energy gap is developed 
which separates off the populated levels. Dots and crosses 
correspond, respectively, to the first and second component. 



in the lower line sum) ones represent the interaction en- 
ergy between atoms within pairs whereas the last term 
describes the pairs interaction (i.e., the interaction be- 
tween distinguishable atoms from different pairs). 

Next, we minimize the expression (|I5ll with respect to 
the real amplitudes w„, Vn under extra condition that the 
total number of pairs equals the number of atoms in each 
component J2n ~ ^^^^ ^^'^^ convenient to 

introduce new variables 0„ in the following way: u„ — 
sinO„ and w„ = cos0„. Now, the energy of the ground 
state Ebcs is a function of angles 0„ and its conditional 
minimum can be found by using the Lagrange multiplier 
technique. One introduces the many-variable function 
Fbcs = Ebcs -AG where G = J2n cos^ 0„ - TV and 
looks for the unconditional minimum of function Fbcs- 
The necessary condition for that turns, together with the 
normalization requirement for the total number of pairs, 
to the following set of equations 



^ ^ ^^n, mm sin 20^^i 2 tan 20^^ ^ ^ ^^m,nm cos O7; 
= {ell'' + + Vnn^nn " A) tan 2e„ 

n = 0,l,2,... 
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FIG. 8: Single-particle Hartree-Fock energies for g — +10.0 
and 3 + 3 atoms system (the upper frame). No energy gap 
is developed. The lower frame shows the same spectrum but 
in the case of 10 + 10 atoms system in the excited state with 
antiferromagnetic phase (as in the lowest part in Fig. 
Surprisingly, the energy gap is visible again although there 
are some states within it. 



tion pattern is visible which means that one component 
serves the other as a periodic potential. This resembles 
somevi^hat the situation when the atoms are confined in 
the optical lattice. For the interaction strength g > —3 
one has Vn ~ l,u„ ~ for n < 3 (and simultaneously 
Vn ~ 0,Un ~ ^ for n > 3) and there is no difference 
between the Hartree-Fock and BCS densities. 
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FIG. 9: BCS energy of 3 + 3 atoms system as a function of 
number of states in the basis of pair states for g = —10. In 
this case Ehf = —23.4 whereas Ebcs ~ —25.7 osc. units. 
The inset shows corresponding occupations in cases when the 
number of basis states equals to 12 (dotted line), 21 (dashed 
line), and 105 (solid line). 



^cos2e™=iV (17) 

rn 

for variables 0„ and the Lagrange multiplier A. We solve 
numerically this set of equations and then check whether 
what we found is a minimum or not. We just randomly 
disturb obtained amplitudes Un,v„ and verify that the 
energy of the system according to expression (I15|l always 
gets higher. 

There is an important issue of choosing the appropri- 
ate two-particle basis functions for pair states. It turns 
out that it is good to build the pair wave functions us- 
ing single-particle Hartree-Fock orbitals. In such a case 
the ground state of the system involves relatively small 
number of pair states. In Fig. Elwe plot the energy of 
the system of 3 + 3 atoms as a function of the coupling 
constant g, calculated within both the Hartree-Fock and 
BCS approaches. It is seen that for strong enough at- 
traction {g < —3) the ground state of the system be- 
comes of BCS type. Simultaneously, the energy gap in 
a single-particle Hartree-Fock spectrum is being devel- 
oped (see Fig. |2l . It is favorable for the system to lower 
its energy by forming pairs of atoms of different com- 
ponents. Fig. |5| shows that, in fact, 50 states basis is 
already good enough. The densities for both fractions 
are plotted in the inset in Fig. El for g ~ —10. In BCS 
case the broad background is developed, not present for 
the Hartree-Fock densities. In both cases the fragmenta- 



According to definition © one can calculate the one- 
particle density matrices for the BCS ground state given 
by in 

n 

P?\yuyi) ^Y.< ¥'l'^*(yi)^l'nyi') (i8) 

n 

Assuming the energetically lowest N pair states are fully 
occupied, i.e. w„<Ar = l,M„<Ar = and Vn>N = 
0, Un>N = 1, the above formulas match the ones derived 
within the Hartree-Fock approach. The same is true also 
for the two-particle density matrix which is given by 

P2^^^^'(xi,yi;x{,y{) =^ (UnWmWnWm + W^W™) X 
n.ni 

V.W*(xi)^,l^)*(yi)^W(xO¥'^,?Hyi') (19) 

In Fig. ^1 we have shown the diagonal part of two- 
particle density matrix for the BCS ground state of 3 + 3 
atoms system for the interaction strength g = —10. For 
this value of g the pairing phenomenon is already well 
established and leads to 10% decrease in energy in com- 
parison with the Hartree-Fock case. The characteristic 
crests along 'x' and 'y' axis are developed as a result of 
presence of higher n (bigger size understood as the mean 
value of |x — y|) pair states in the BCS ground state (see 
also Fig. EJ. No such crests are visible for the Hartree- 
Fock ,, ground state". 
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FIG. 10: Diagonal part ol two-particle density matrix for the 
BCS ground state of 3 + 3 atoms system, i.e., joined proba- 
bility of detecting two distinguishable atoms. The interaction 
strength g = —10. All quantities are given in oscillatory units. 
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FIG. 12; Single-particle Hartree-Fock energies for various 
strengths of attraction between atoms of different components 
for a system of 50 + 50 atoms. For strong enough attrac- 
tion the energy gap is developed (upper frame) and the BCS 
ground state turns out to be reachable. When the interaction 
strength is increasing further the energy gap is getting bigger 
(lower frame). 
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FIG. 11: Correlation function for the BCS ground state of 
3-1-3 atoms system. The interaction strength g = —10. Only 
this part of 'xy' plain is considered where the two-particle 
density matrix essentially differs from zero. All quantities are 
given in oscillatory units. 



In Fig. 111! we plot the correlation function g2, defined 
generally as 



52(x,y) 



P2^^*^^ (x, y; X, y) - p\'' (x; x)pi^^ (y; y) 



.(2) 



(2), 



(20) 



in the case of one-dimensional 3 + 3 atoms system un- 
der the same parameters as in the case of Fig. ^1 It is 
true that the correlation function differs from zero mainly 
along the directions y = x and y — ~x. Whereas the for- 
mer is not surprising because densities of both fractions 
look the same (see inset in Fig. O, the latter means that 



the strong correlations are developed for distinguishable 
atoms being at positions (2,4) and (—4,-2) in oscilla- 
tory units. Since it corresponds to the size of pair states 
contributing essentially to the ground state of 3-1-3 atoms 
system (approximately first 50 pair states - see Fig. ^ it 
means that the presence oi y = —x correlations supports 
the explanation of properties of the ground state with 
the help of pairing phenomenon. 

Increasing the number of atoms in both fractions we 
find further qualitative changes in the properties of the 
ground state in comparison with 3-1-3 case. When the 
coupling constant g becomes smaller (i.e., for stronger 
attraction) the energy gap in the Hartree-Fock spectrum 
gets bigger (Fig. [T^ just like in the case of 3 -|- 3 system. 
However, some new features appear when one looks at the 
size of the pairs (Fig. I13|l as well as at their energies (Fig. 
I14II for various values of the coupling constant g. For 
strong enough attraction the size of the pairs is getting 
frozen revealing plateau and the gap is developed in the 
energy spectrum. The energy gap in the pairs spectrum is 
perhaps the signature that the system enters a new phase 
characterized by the appearance of clusters of particles 
in its ground state, each built of four atoms. 

It remains to be checked what is the range of valid- 
ity of parameters used in our calculations. To this end, 
we show the following simple estimation. Binary colli- 
sions dominate provided the mean interparticle separa- 
tion is much bigger than the s-wave scattering length, 
i.e., ri3_D|a3_Dp <C 1. According to Olshanii's result |2(]| |. 
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FIG. 13: Size of the pair states calculated as the mean value 
of [x — yl for the 50 + 50 atoms system. Note that for strong 
enough attraction the size is getting frozen and characteristic 
plateau is developed. 
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FIG. 14: The energy of the pair states for the system of 50+50 
atoms. For strong enough attraction the gap in the pairs 
spectrum is developed. 



for elongated systems the efFective one-dimensional in- 
teraction strength can be introduced via the relation 



giD 



930 



1~C 



0-30 



(21) 



where g^D = iTrh^a^u/m, C — 1.4603, and a±_ = 
{h/mLOi_Y/'^ determines the radial length scale. As- 
suming the tight confinement one can introduce one- 
dimensional density in such a way that 



The range of validity of our model can be then expressed 
as 



7r|5ii5Cai/2 + 4|3 



«1, 



(23) 



niD 



(22) 



where a = uJz/(^j_ and quantities nio and giD are already 
taken in units related to one-dimensional system. For 
example, for \giD\ = 10 the condition H23() leads to the 
relation a <C 0.075 (or oj±/uJz ^ 13) which, first, acts 
in the same direction as the assumptions leading to the 
formula H21() and, second, can be fulfilled experimentally 
although for larger \giD\ the radial squeezing is getting 
stronger and stronger. 



VI. CONCLUSIONS 

We have analyzed the system of two-component Fermi 
gas confined in a magnetic trap at zero temperature un- 
der assumption that the only interparticle interaction left 
is the repulsion or attraction between atoms of different 
components. We employ the microscopic description of 
the system by using the explicit atomic wave functions. 
The structure of the ground state strongly depends on 
the sign of the coupling constant. When distinguishable 
atoms repel each other strongly enough both fractions 
separate with that which is confined stronger filling the 
center of the trap. In the opposite case, i.e., when the 
atoms strongly attract each other, the ground state dra- 
matically changes its character. The BCS pairs appear 
which is proved by showing that the energy of the sys- 
tem of pairs is lower than the corresponding Hartree- 
Fock energy and by analyzing the two-particle correla- 
tion function. For even stronger attraction the system 
enters a new phase where presumably clusters of parti- 
cles are formed. We are working currently on extension of 
these calculations to the three-dimensional case as well 
as on inclusion of the optical lattice and finite tempera- 
tures. 
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